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A pulse of light, injected into a weakly disordered dielectric medium, typically, will leave its 
' initial location in a short time, by diffusion. However, due to some rare configurations of disorder, 

, there is a possibility of formation of high quality resonators which can trap light for a long time. 

' We present a rather detailed, quantitative study of such random resonators and of the "almost 

^ ' localized" states that they can support. After presenting a brief review of the earlier work on the 

^ , subject, we concentrate on a detailed computation of the "prefactor": knowledge of the latter is 

crucial for varifying the viability of the random rasonators and their areal density. Both short 
range disorder (white noise) and correlated disorder are studied, and the important effect of the 
■ correlation radius, Rc, on the probability of formation of resonators with a given quality factor Q is 

discussed. The random resonators are "self- formed", in the sense that no sharp features (like Mie 
scatterers or other "resonant entities") are introduced: our model is a featureless dielectric medium 
(— I ■ with fluctuating dielectric constant. We point out the relevance of the random resonators to the 

' recently discovered phenomenon of coherent "random" lasing and review the existing work on that 

subject. We emphasize, however, that the random resonators exist already in the passive medium: 
gain is only needed to "make them visible". 
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I. INTRODUCTION 



I There exists a formal analogy between the Schrodinger equation describing electron motion in a random potential, 
(~| : and the scalar wave equation for light propagation in a medium with fluctuating dielectric constant. Consider, for 
O ' concreteness, a two-dimensional geometry (quantum well for an electron and a thin film for light). Then both equations 
can be presented in the form 



Ap^ + - Uip)] = 0, (1) 



where p is the in-plane coordinate. For an electron with mass m and energy E, moving in a random potential V{p), 
the parameters fc^ and U (p) are 

^ ; fc^ = 2mE , U{p) = 2mV{p). (2) 

For a light wave with frequency w, traveling in a medium with dielectric constant e, corresponding expressions for k 
, and U (p) take the form 

fc^-f-)^ Uip)^-Seip)n\ (3) 



c / \ c 



where Se{p) is the fluctuating part of the dielectric constant. 

The analogy between electrons and light is incomplete. For one thing, the "scattering potential" for light depends 
Q , on frequency and it vanishes in the limit of zero frequency. Moreover, in a dielectric medium (positive e) must be 
O ' positive which implies that there can be no bound states for light, as long as Se{p) is assumed to vanish outside some 
^ finite region in space. Therefore the simple notion of a binding potential well does not exist for light: for instance, 
a dielectric sphere embedded in a uniform medium cannot bind photons, regardless of whether its dielectric constant 
is larger or smaller than that of the surrounding medium. Thus, unlike electrons, photons will always escape from a 
dielectric medium into the surrounding air. However, under appropriate conditions, they can be trapped within the 
sample for a long time. The main purpose of this paper is to discuss trapping of light in a weakly disordered dielectric 
film. 

The random term in Eq. (|l|) is zero on average; its statistical properties are described by the r.m.s. value, Uq, and 
the correlator K{pi — P2) 

{U{p)) = , {U{p^)U{p2)) = UlK{p^ - P2) , K{0) = 1. (4) 

It is known Q that in two dimensions the nature of solutions of Eq. (|l]) is governed by the dimensionless conductance 
kl, where I is the transport mean free path. Using the golden rule, the product kl can be expressed through the 
correlator, K, as follows 
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(kl)-^ = J dq d(j) q^S{q^ + 2fcgcos0) J d^p K{p) exp(igp) = 

/'7r/2 />oo 

/ da sin^ a / dp p K{p) Jo(2/cpsina), (5) 
Jo Jo 
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where Jq is the Bessel function of zero order. The above integral can be evaluated analytically if we choose a gaussian 
form for the correlator K{p) = exp(— p^/i?^). Substituting this form into (^), we obtain 

where the dimensionless function F is defined as 

F(x)=e-^[/o(x)-/i(a;)]. (7) 

Here Iq and Ii are the modified Bessel functions of zero and first order, respectively. 

As a simple example of a realistic two-dimensional disorder, consider a system of disks with random positions of 
the centers and with fluctuating radii described by a normalized distribution function, Within a disk, we have 

U{p) = Ud, whereas outside the disk U{p) ~ 0. If the filling fraction, /, is low, / ^ 1, the positions of the centers 
of the disks are uncorrelated. Then K{p) is determined by the overlap between two circles of the same radius with 
centers shifted by p. A straightforward calculation yields 



(C/(pi)C/(pi +p))^^JdR e(R 



^ V2i?y 2R V \2r) j ' 

where 9{x) is the step-function. The actual distribution function, is governed by various technological factors. 

However, it is reasonable to assume that small values of R are strongly unlikely, and that ^{R) falls off abruptly 
at large R. Consider as an example the distribution <i>(x) = c'^x'* exp(— cx^), where x — R/R {R is the average 
radius) and c = 64/97r. The distribution is designed so as to yield the value 50% for the relative spread in R, i.e., 
{\R — R\)/R = 0.5. The result of calculation of the integral (g) using this distribution is shown in Fig. 1 together 

with its gaussian fit, which yields C/q = 0.86C/J/ and K{p) = exp (-3Ap'^/R^y From this fit we conclude that the 

gaussian form of the correlator, K[p), corresponds to a quite generic distribution of R. 
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Equation (^) provides the explanation why the criterion for strong localization, kl < 1, can be easily satisfied for 
electrons, but very hard to achieve in the case of photons. Indeed, for electrons Eq. (|^) can be presented as 

ikl)-'-{^)Vo'RlF{mERl), (9) 

while for photons Eq. takes the form 

(H,-..(iL)(^)'A»«?f(f^), a„, 

where Vq and A are the r.m.s. fluctuations of the potential and dielectric constant, respectively. It is seen from Eq. 
(|) that, since the function F is always smaller than one, the electrons are strongly localized in the low-energy domain 

E < Ec- The value of Ec is different for the short-range, Rc ^ (mVb)~^^^, and for the smooth, Rc » (toVo)"^^^, 

potentials. For short-range potential we can use the asymptotics ~ 1, which leads us to the estimate 

Ec ~ rnVgRl- Hence, mEcR^ ~ (mVb-Rc)^ ^ 1' that the argument of F is indeed small at E = Ec- In the case of 
the smooth potential, the condition mVoRl ^ 1 suggests that the semiclassical description applies, so that Ec ^ Vc 
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Calculation of the mean free path based on the golden rule is inadequate in this case. To see this, note that at E' ~ Vq, 
the argument of F in Eq. is large, so that we can use the asymptotics oc Then Eq. yields 

I ~ Rc for E ~ Vq. Thus, for smaller E, namely E < Vq, we have I < Rc- The latter relation indicates the failure of 
the perturbation theory for E <Vo. 

Now let us perform the similar analysis for photons. Using the small-x and large- a; asymptotics of the function F, 
we obtain from Eq. ( p^ ) 

= nikoRcAr ' = 7ri/2A2 ' 

where feg — uj/c. Eq. ( |ll|) suggests that when the r.m.s. fluctuation A is weak, A ^ e, we have kl ^ I both in the 
low-frequency (k^Rc «C 1) and in the high-frequency {k^Rc ^ 1) domains. This peculiar result is due to the already 
mentioned frequency dependence of the "optical potential". 

The above analysis, however, does not rule out the possibility of light localization in the strongly scattering media. 
In fact, first experimental indications of localization effects for microwaves were reported more than a decade ago [^^. 
In these experiments localization was inferred from the measurements of various transmission characteristics of the 
microwaves through the tube filled with a random mixture of aluminum and Teflon spheres. For optical frequencies 
l^,^!^, the strongly scattering medium used in first experiments, aimed at light localization, was a semiconductor 
{GaAs) powder. Measurements of the transmission vs. the sample thickness were complimented with measurements 
of the coherent backscattering (CBS) cone |@,||. In the latter measurements localization manifests itself through the 
rounding of the CBS cone by limiting the maximal length of coherent path 

The early reports of the observation of wave localization both for microwaves and for light [2-4] were inconclusive 
because of the possibility that the results were affected by absorption. To get rid of this ambiguity, in the later 
experiments the CBS mesurements [|o| were performed on the macroporous GaP networks, which scatter light stronger 
than a powder . This allowed the authors to rule out the absorption or the finite sample size as a source of rounding 
of the CBS cone. For microwaves [11-13], the recent progress in detecting localization is due to a novel approach 
to the analysis of the transmission data based on analysis of the relative size of the transmission fluctuations. This 
approach permits one to detect localization even in the presence of absorption. 

In the weakly scattering active medium the propagation of light remains diffusive. However, the interplay of the 
diffusion and the gain-induced amplification can be very nontrivial. Namely, this interplay can give rise to the 
incoherent random lasing, predicted by V. I. Letokhov [Q. As it was pointed out in Ref. jlj], there is a close 
analogy between multiplication of neutrons in course of the chain reaction and photons in the amplifying disordered 
medium (photonic bomb). Upon first experimental observation of incoherent random lasing [p^ , it was subsequently 
reproduced for various realizations of the gain media and different types of disorder. Comparison of theoretical 
[16-21] and experimental results [15,22-42] has confirmed that the diffusion theory, which neglects the interference 
effects, is quite sufficient for the description of incoherent lasing. Except for studies on powder grains of laser 
crystal materials [37-40] and 7r-conjugated polymer films |4^,^^, the majority of experiments [15,22-36] have used 
dye solutions as amplifying media. Colloidal particles suspended in a solution served as random scatterers. These 
scatterers are responsible for nonresonant feedback required for incoherent lasing. The essence of this feedback is that 
the light amplification length, la (in the absence of disorder), is significantly shortened [to ^ {lla)^^'^] in the disordered 
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medium, when the light propagation is diffusive. The threshold condition for incoherent lasing corresponds to the 
gain magnitude at which (/?a)^^^ becomes of the order of the sample size. 

In contrast to incoherent lasing, the recent discovery of coherent random lasing [43-46] adds a new dimension to the 
physics of light propagation in disordered media. Coherent random lasing emerges as the degree of disorder increases, 
so that the mean free path, I, becomes progressively smaller. The fact that the light, emitted from a disordered 
sample, is truly coherent, which does not necessarily follow p^ , ^8|] from drastic narrowing of the emission spectrum, 
observed in Refs. [43-46], was later convincingly demonstrated in photon statistics experiments [^,Q. In the absence 
of mirrors, it is evident that, in order to support the coherent lasing, the disordered medium itself should assume 
their role. The latter is feasible only due to the interference effects, that are not captured within the diffusion picture. 
More quantitatively, in order for the random medium to play the role of a Fabry-Perot resonator, it is necessary that 
certain eigenfunctions of Eq. (|l|) were either completely localized or almost localized. Almost localized solution can 
be, roughly, envisaged as a very high local maximum of the extended eigenfunction ip{p) of Eq. (|l]). If this maximum 
is viewed as a core of ip{p), then the delocalized tail (see Fig. 2) can be viewed as a source of leakage. In other words, 
the core itself, being not an exact eigenfunction, can be viewed as a solution of Eq. corresponding to a complex 
eigenvalue Imfc^ ^ 0. Then the weak leakage translates into a small value of the imaginary part of . Inverse of this 
imaginary part determines the lifetime of the core. The higher is the local maximum of ^{p), the longer is the lifetime. 
Making link to the Fabry-Perot resonator, the ratio fc^/Imfc^ can be identified with a quality factor, Q. Thus, from 
the perspective of the coherent random lasing, the right question to be asked is: how high are the attainable quality 
factors of the almost localized solutions of Eq. (0) at a given prameters of disorder, i.e., magnitude, A, and correlation 
radius. Re- 



(a) 




FIG. 2. (a) Spatial distribution of the wave function of the anomalously localized solution, 'ip{p), of Eq. (^) is shown 
schematically, (b) Spatial distribution of the dielectric constant, e(p), corresponding to the trap, responsible for the solution 
e(p) = e outside the blank region. Only the lower half of the trap is shown. 

We address this question in the present paper. However, prior to discussing the almost localized solutions of Eq. 
(|l|), a subtle point must be clarified. Namely, Eq. (|l|) with the "potential" being real, describes propagation of light in 
a passive medium. It might be argued that, in the presence of gain, which is required for lasing, the spatial structure 
of the solutions of Eq. (|l|) undergoes a drastic change, so that the almost localized modes of a passive random 
medium and actual lasing modes have little in common. In fact, it is well known both from theory [ pT[ ] and from 
CBS experiments that the diffusive trajectories of light within a random medium get elongated in the presence of 
the gain. With regard to coherent random lasing, it was initially claimed [5^ that the gain facilitates localization of 
light. This claim was even supported by the numerical simulations |5^. However, in the later theoretical [54-56] and 
experimental ]^7| papers it was explicitly stated that, similarly to the conventional lasers with Fabry-Perot resonators, 
the gain only reveals high-Q solutions of Eq. (|l|) existing in the passive medium. Thus, in the present paper, we will 
focus exclusively on the passive disordered media. The question about the likelihood of formation of the disorder- 
induced resonators is central to the understanding of the coherent random lasing. This question is at the core of the 
ongoing in-depth experimental studies [57-65]. Except for Ref. [^, theoretical papers on random lasing [56,67-71] 
do not address this question. 

Let us finally mention that a similar question of trapping electrons, for a long time, in a weakly disordered conductor 
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was addressed long ago in the context of electronic transport [72-75]. However, the "almost localized" states discussed 
in 1^,^,0 and, further, in this paper are quite different from the "prelocahzed" states studied in [72-75]. Wave 
functions of the "almost localized" states are confined primarily to a small ring (a high-Q resonator). Moreover, 
these states are extremely sensitive to the value of the correlation radius, Rc, of the disordered potential. We shall 
see below that, for a given parameter fc^ > 1, the likelihood of high-Q "almost localized" modes is very low for the 
short-range disorder but sharply increases with Rc- Such features are absent for the "prelocahzed" states [72-75], with 
their comparatively large spatial extent. 



II. LIKELIHOOD OF RANDOM RESONATORS IN A DISORDERED FILM 



A. Intuitive Scenario of Random Cavities Based on Recurrent Scattering 

The first experimental work on coherent random lasing by Cao et. al. [^3| was carried out on thin (with a width 
~ 27r/fco) zinc oxide (ZnO) polycrystalline films. Laser action manifested itself through a drastic narrowing of the 
emission spectrum when the optical excitation power exceeded a certain threshold. The authors of Ref. |^3| realized 
at the time that coherent lasing requires a resonator (cavity). They conjectured that, in the absence of traditional 
well-defined resonator (as in a semiconductor laser), the cavities in polycrystalline films are "self- formed" due to strong 
optical scattering. For a microscopic scenario of the formation of such cavities, they alluded to the remark made in 



Ref. 1 18 1 that closed-loop paths of light can serve as "random ring cavities". The importance of these closed-loop 
paths (recurrent scattering events) was pointed out earlier ||78| , when they were invoked for the explanation of the 
magnitude of the CBS albedo. With regard to CBS, the effect of recurrent scattering events, e.g., the events in which 
the first scattering (of the incident wave) and the last scattering (of the outgoing wave) are provided by the same 
scatterer, is that these events do not contribute to the CBS, thus reducing the albedo in the backscatering direction. 

The fact that recurrent events show up in CBS certainly does not allow to automatically identify closed-loop paths 
with resonator cavities. Therefore, the feasibility of the scenario of random cavities, adopted in Refs. |^,^,^,^ for 
interpretation of experimental results, was later put in question |79| . The arguments against this scenario were the 
following. Since in each scattering act most of the energy gets scattered out of the loop, an unrealistically high gain 
would be required to achieve the lasing threshold condition for such a loop. Also the loops of scatterers are likely to 
generate a broad frequency spectrum rather than isolated resonances. 

Certainly the picture of random cavities representing a certain spatial arrangement of isolated scatterers is too 
naive. This, however, does not rule out the entire concept of disorder-induced resonators. Although sparse, the 
disorder configurations that trap the light for long enough time can be found in a sample of a large enough size, and a 
single such configuration is already sufficient for lasing to occur. Therefore, under the condition kl ^ 1, which implies 
that overall scattering is weak, the conclusion about the relevance of random cavities can be drawn only upon the 
calculation of their likelihood. This calculation is decribed in the rest of this Section first on qualitative, and then on 
quantitative levels. 

Making link to the discussion in the Introduction: by pursuing the random cavity scenario Cao et. al. have 
intuitively arrived to the concept of prelocahzed states, that was introduced in transport more than a decade ago. 
The fundamental difference between the prelocahzed states and recurrent events is that recurrent events emerge in 
the higher order [in parameter kl~^] of the perturbation theory (they correspond to the specific type of diagrams called 
Hikami-boxes [pO[), whereas the formation of prelocahzed states is a genuinely nonperturhative effect. Already in the 
first analytical approach to the problem of prelocahzed states [ ^ , it was demonstrated that, in order to capture an 
anomalously slow tail of the conductance relaxation (which is due to trapping), all the orders of the perturbation 
theory must be taken into account. 



B. Optimal Fluctuation Approach to the Problem 



1. Qualitative Discussion. 

Similarly to the treatment in p^ , |5^ , |59| we restrict our consideration to the two-dimensional case (a disordered 
film). Regarding the geometry of a random resonator, we adopt the idea proposed by Karpov for trapping acoustic 
waves |7^ and electrons [^^ in three dimensions. 

Suppose that within a certain stripe the effective in-plane dielectric constant of a film is enhanced by some small 
value ei <C e. Then such a stripe can play a role of a waveguide^ i.e., it can capture a transverse mode, as it is 
illustrated in Fig. 3. There is no threshold for such a waveguiding, which means that the transverse mode will be 
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captured even if the width of the stripe is small. Now, in order to form a resonator, one has to roll the stripe into a 
ring. Upon this procedure, the mode propagating along the waveguide transforms into a whispering-gallery mode of a 
ring. An immediate consequence of the curving of the waveguide is emergence of the evanescent leakage - the optical 
analog of the under-the-barrier tunneling in quantum mechanics (see Fig. 2). This leakage is responsible for a finite 
lifetime of the whispering-gallery mode. Thus we have specified the structure of the weakly decaying solutions of Eq. 
(|l|), discussed in the Introduction. Namely, the mode of the waveguide plays the role of the core, while delocalized 
tail (see Fig. 2) reflects the evanescent leakage. Due to the azimuthal symmetry, the modes of the resonator are 
characterized by the angular momentum, m. Denote by Nm(fc^, Q) the areal density of resonators with quality factor 
Q in the film with a transport mean free path I. Obviously, in the diffusive regime, kl > 1, the density N,„(fc^, Q) is 
exponentially small for Q 3> 1. In this domain N,„(fc^, Q) can be presented as 

N„,(M,g)cxe-^"('^'''3) . (12) 

Let us first give a qualitative estimate for Smi which reveals its sensitivity to the strength and the range of the 
disorder (A^ and Rc). Since m is the number of wave lengths along the ring, its radius is po = m/k. The ring 
waveguide can support a weakly decaying mode only if its width w satisfies the condition w(ei/e)^^^ > k^^. A 
straightforward estimate for the decay time due to evanescent leakage, i.e., the quality factor Q of the waveguide 
results in hiQ ^ fcpo(ei/e)^^^. Since the number kpo = m is large, a relatively small fluctuation of the dielectric 
constant within the area 2ttpqw of the ring can produce a large value of Q. 




FIG. 3. Rationale for the structure of the resonator. Upon wrapping a stripe with enhanced dielectric constant into a ring, 
a waveguided mode transforms into a whispering-gallery mode. 

The probability W for creating the required fluctuation strongly depends on Rc- For a short range disorder 
(k^Rc ^ 1) fluctuations of order ei should occur independently in a large number, N ^ pqw/R^, of spots within 
the ring, so that the probability W ~ exp (— A^e^/A^). In the other extreme of strongly correlated disorder, when 
Rc ^ w, the number of independent spots is much smaller, N' ~ Po/Rc (the number of squares with a size Rc needed 
to cover the ring). Correspondingly, the probability W ~ exp (— 7V'e^/A^) is much larger than for the short range 
case. Finally, using the relation between ei and Q, the probability W = exp(— S*™) can be rewritten in terms of Q, 
thus, yielding an estimate for Sm- For the short range case {k^Rc <^ 1) we obtain Sm ~ kllnQ, where the mean free 
path I is proportional to (i?cA)~^ [see Eq. (pd])]. In the opposite limit of a smooth disorder Eq. ( pd| ) yields for the 
transport mean free path / - Rc/A"^. Then we have 5*^ ~ TV'ef/A^ ~ Z(lnQ)4/3/(fci?2^V3)^ Thus, for given kl and 
Q, the density of resonators is the higher the smoother is the disorder. This conclusion is central to our study and 
will be addressed below in more detail. 



2. Quantitative Results from the Optimal Fluctuation Approach. 

The above program can be carried out analytically [|66| with the use of the optimal fluctuation approach |^ , ^ . 
This approach is based on the idea that, when the exponent, Sm, in Nm(fc?, Q) is large, then the major contribution 
to Nm(fcZ,(5) comes from a certain specific disorder realization. In application to random resonators, the optimal 
fluctuation procedure reduces to finding the most probable fiuctuation of the dielectric constant which is able to trap 
the light for a long time ^ u^^Q. Assuming that the fluctuation is azimuthally symmetric (see Fig. 2), the shape of 
the optimal fluctuation can be found explicitly [ ]66[ yielding the following expression for exponent Sm 
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where ei = e(31n(5/2TO)^/^. The analytical expression for the function $(71) is the following 

^, , 33/2 (5 + V9 + 16^2)5/2 

= — — ^ — . (14) 

^ ^ 26 (3 + x/9 + 16^2)3/2 ^ ^ 

Recall now, that we are interested in the density of random resonators at a given value of kl. The transition from A 
to / is accomplished by using Eq. (11). For the short range case, when Rc and $(ii) 1, we obtain 

^\ 1/2 



SmikoRc <^ 1) ^ 2 [ — ] kllnQ . (15) 



To trace the change of Sm with increasing Rc it is convenient, after using Eq. (|1 1\) , to present Eq. ( |13D in the form 

SmikoRc > 1) _ HeY^koRc) 



SmikoRc « 1) 7ri/2(ei/2fcoi?^)3 



(16) 



It is seen from Eq. ( p^ ) that S'„i falls off rapidly with increasing Rc- In the domain koRc > 1, but e^^koRc < 1 we 
have $ « 1, so that Sm cx (fcoi?c) For larger Rc we have $(u) c>c u. In this domain Sm decreases slower with Rc'. 

Sm oc {koRc)^ 

5™(^oi?.»l) = ^^;^,^^. (17) 

Asymptotic expressions (p^ ) and ([l7| ) agree with the results of the qualitative derivation, with all numerical factors 
now being determined. We emphasize that Eqs. ( p^ ) and ( p^ apply for a given kl value, so that the decrease of Sm 
with Rc leaves the backscattering cone unchanged. 



3. Estimates. 

Equation (|l5|) qua ntifies the effectiveness of trapping of light in a random medium with point-like scatterers. It 
follows from Eq. (ttSh that the likelihood of high-Q cavity is really small. Indeed, even for rather strong disorder, 
kl = 5, the exponent, 5*^, in the probability of having a cavity with a quality factor Q = 50 is close to Sm = 120. We 
emphasize that in two dimensional case under consideration, this exponent does not depend on m and, thus, on the 
cavity radius po = ni/t^/'^kQ. More accurate calculation [Q, taking into account the corrections to Eq. (|l5|), indicates 

1/2 

that Sm as a function of m has a minimum at m ^ (kl In Q) . 

To estimate the degree to which finite size of scatterers (^^ Rc) improves the situation, we choose k^Rc ~ 2, which 
already corresponds to the limit fco-Rc S> 1 in Eq. (pi]), but still allows to set <i> = 1. Then for Q — 50, kl — 5 we 
obtain 5*^ ~ 1.1, suggesting that the resonators with this Q are quite frequent. In the latter estimate we have set 
e = 4. 



C. Frequently Asked Questions 

1. Why Rings? 

The answer to this question is illustrated in Fig. 4. The distinguishing property of a ring is that the local curvature 
radius is the same at each point. Upon any deviation from the ring geometry, the curvature in a certain region of 
the fluctuation would be higher than in all other regions. Since the evanescent losses are governed by this curvature, 
the quality factor of the resonator would be determined exclusively by this region (see Fig. 4) , so that the remaining 
low-curvature part would be "unnecessary" , in the sense, that a ring with a radius corresponding to the maximal 
curvature would have the same quality factor as a square in Fig. 4 but significantly higher probability of formation. 
It is also quite obvious that, for the purpose of supporting a wave-guided mode of the whispering-gallery type, a ring 
is much superior to a disk of the same radius: indeed, the internal area of the disk remains unused in the guiding 
process, whereas a heavy penalty in terms of probability is paid in creating this area. 
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evanescent leakage 

FIG. 4. This drawing illustrates the optimal character of the ring-shape resonator. 



2. Why Smooth Disorder Facilitates Trapping? 



At the qualitative level, the enhancement of the probability of formation of the cavity with increasing Rc can be 
understood for a toy model of the disorder, illustrated in Fig. 5. Suppose that all the disks, that model the scatterers, 
are identical. Then Rc scales with the radius of the disk, R. Since the disks cannot interpenetrate, the ring-shaped 
cavity corresponds to their arrangement in the form of a necklace. The probability of formation of such a cavity 
can be estimated as follows. Suppose that a sector, (50, is "allocated" for a single disk. The probability to find a 
disk within this sector, at the distance po from the center , is ^ n(po^0)^, where n is the concentration of the disks. 



Thus, the probability of formation of the necklace is exp 



■In 



where I?- is the number of sectors. It 



is obvious that if a necklace is "loose" , the quality factor of the corresponding cavity would be low. In order for Q to 
be high, neighboring disks must almost touch each other. This implies that d(f> ~ 2R/ pq. Then the above estimate 



for probability takes the form exp 



-^ln(i) 



where / = mrR^ is the filling fraction. This probability increases 

exponentially with R, i.e., with Rc, reflecting the fact that, for a given po, the number of disks to be arranged is smaller 
when R is larger. The above estimate was based on the assumption that the positions of the disks are uncorrelated, 
i.e., / ^ 1 (in contrast to |Q where / — 0.4). We have used the model of hard disks as an easiest illustration of 
the role of Rc- Obviously, Eq. (|l^) does not apply to this model, since it was derived under the assumption that the 
statistics of the fluctuations is gaussian. 



Po 




FIG. 5. A schematic illustration explaining why larger correlation radius for a fixed filling fraction facilitates trapping. A 
sector, shown with dashed lines, illustrates the tolerance in the arrangement of disks into a necklace. 
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3. "Vulnerability" of the Ring-Shaped Cavities. 

The value Sm given by Eq. (|l3|), which was derived within the optimal fluctuation approach, is the exponent in 
the probability of formation of an ideal ring. Obviously, any actual disorder realization is not ideal, in the sense, that 
actual distribution of dielectric constant differs from the optimal. For the same reason, the probability of formation 
of ideal necklace of the type shown in Fig. 5 is zero. In order for the probability to be finite, we should allow a certain 
tolerance in the positions of the centers of the disks, as illustrated in Fig. 5. In the conventional applications of the 
optimal fluctuation approach deviations from the optimal distribution do not affect the value of the exponent, 
Sm- However, in application to the random cavities, we have searched for the fluctuation which is optimal for a given 
trapping time, lu^^Q. In this particular application, a "normal" gaussian deviation from the optimal geometry can 
have a catastrophic effect on trapping by scattering the light wave out of the whispering-gallery trajectory. This 
scattering is discussed below. 

Scattering within the plane. Two-dimensional picture adopted throughout this paper, implies that electromagnetic 
field is confined within a thin film in the z-direction. This confinement results from the fact that the average dielectric 
constant of the film is higher than in the adjacent regions. Then the filed distribution, £o{z), along z corresponds to 
a transverse waveguide mode. For a given frequency, uj, the almost localized state on a ring can be destroyed due 
to the scattering into states with the same distribution of the field in the z-direction, which propagate freely along 
the film. More precisely, the almost localized state with a given angular momentum to, which is protected from the 
outside world by the centrifugal barrier, can be scattered out to the continua of states with smaller to's, for which 
there is no barrier. It is essential to estimate the lifetime, r, with respect to these scattering processes and to verify 
that it is feasible to have r larger than the prescribed trapping time, lu~^Q, so that the almost localized state is 
not destroyed. A rigorous treatment of this "scattering out" effect is quite involved and is done in Section III, for a 
Gaussian potential. 

The effect can be also illustrated with the model of randomly positioned hard disks (Fig. 5), although the disorder 
in this model is non-gaussian. It is seen from Fig. 5, that spacings between the rings, which are due to tolerance, 
open a channel for the light escape, that is different from evanescent leakage. A typical lifetime with respect to such 
an escape is quite short, i.e., even a small tolerance, which affects weakly the exponent in the probability of the cavity 
formation, seems to be detrimental for trapping. At this point we emphasize that, in calculation of the scattering rate 
out of the whispering-gallery trajectory, the disks constituting the necklace must be considered as a single entity. As 
a result, for a given configuration of the disks, the rate of scattering out caused by a single disk must be multiplied 
by the following form-factor 



T = 



27r 



exp {ikpi 



^Jo(/c|p,-p,|), (18) 



\ 2 



where pi is the position of the center of the z-th disk in the necklace. The form-factor, JF, is the sum of (j-^) S> 1 
terms. Out of this number, {kR)~^ (■^) terms (for kR < 1) and terms (for kR > 1), for which k\pi — pj\ < 1, 
are close to unity. The portion of these terms is small. Other terms have random sign. This leads us to the important 
conclusion that, for certain realizations of the necklace in Fig. 5, the form-factor can take anomalously small values. 
For these realizations the quality factor will be still determined by the evanescent leakage. The "phase volume" of 
these realizaions is exponentially small and depends strongly on the model of the disorder. 

Scattering out of the plane. Compared to the previous case, two modifications are in order. Firstly, since the final 
state of the scattered cavity mode is a plane wave with the wave vector pointing in a certain direction within the solid 
angle in, the expression Eq. ( p^ ) for the form-factor should be replaced by 

^^^sin(fco|p.-p,|)^ (19) 
ko\p.-pA 

Secondly, for kR > 1, i.e., when the disorder is smooth, scattering out of the plane that is caused by a single disk, 
requires a large wave vector transfer, ^ k. Thus, the corresponding rate is suppressed as compared to the in-plane 
scattering. 
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III. PREFACTOR 



A. Qualitative discussion 



In Sec. II we calculated the probability of the formation of a trap that captures hght for a long time r = uj~^Q. 
However, the relevant characteristics of traps is their ureal density. Within the optimal fluctuation approach, the 
relation between the probability and the areal density emerges in course of calculation of the prefactor [ p5p^ ] . Namely, 
the combination with dimensionality of the inverse area comes from the so called "zero modes" , which reflect the fact 
that the fluctuation can be shifted as a whole in both x and y directions. A typical shift, making two fluctuations 
independent, is of the order of the extent of the fluctuation in each direction. This suggests that the proportionality 
coefficient between the density of traps and the probability of the formation of the trap is roughly the inverse area 
of the fluctuation |8a] . In our particular case, when the fluctuation is ring-shaped, the estimate for the dimensional 



prefactor is 



where w is the width of the ring (see Sec. 11). The dimensionless part of the prefactor within the 



optimal fluctuation approach reflects the "phase volume" of the fluctuations, which perturb the shape of the optimal 
fluctuation leaving the "energy" fc^ unchanged. 

For the almost localized modes, considered above, the situation with prefactor is qualitatively different from the 
case of the truly localized states, for which the optimal fluctuation approach was devised The specifics of the 

almost localized states is that their "energy", A:^, [see Eq. (|l|)] is degenerate with continuum of the propagating modes. 
As a result, a typical small perturbation of the shape of the ring will not only shift k^, but also cause the coupling 
of the trapped mode to the continuum, or in other words, the additional leakage will emerge due to the fluctuations. 
In order to incorporate this effect into the theory, the density of traps should be multiplied by the probability, P{t), 
that the lifetime with respect to this additional leakage is longer than r |84|. To estimate this probability, we consider 
the additonal leakage for a given disorder realization 



dp %{p)U{p) ^f.{p) 



d{k 



(20) 



where tpQ is the "localized" solution of Eq. (Il|) (without evanescent leakage) , and k^ are the propagating eigen- 



functions and eigenvalues of Eq. 
small matrix elements (V'o|t^|V'Ai) 



(H), respectively. It is seen from Eq. (|2C|) that small additional leakage requires 
Since the functions belong to the continuous spectrum, it might seem that this 



requirement is impossible to meet. This is, however, not the case, since normalization factor of \ipfj.\ is the inverse 



system size. To demonstrate that additional leakage can be small, 
namely, the plane waves, and rewrite Eq. (20) in the form 



it is convenient to use the explicit form of tp^ 



~ = JJ dpidp2U{pi)Sipi,p2)U{p2), 



(21) 



where the kernel S{pi, P2) is defined as 



S{pi,P2) = Vo(Pi)V'o(P2) (^) 



27r 



2k 



^o(Pi)V'o(P2)^o {k\Pi - P2I) , 



(22) 



It is seen from Eq. ( [2^ ) that the kernel S{pi,p2) is exponentially small if one of the points, pi or p2, is located 
outside the region, occupied by the "body" of ■00 (p)- When both pi and p2 are located inside this region, then the 
characteristic spatial scale of the kernel, S, is |pi — P2I ~ Thus, for the sake of our qualitative discussion, we 

can replace S{p\,p2) by A~^6 {k~^ — \pi — P2I), where A ^ |^o(0)|^^ is the area of the fluctuation, and restrict 
integration in Eq. ( |2l| ) to the region of the area A. 

Averaging over disorder configurations in Eq. (|l|) yields the mean value of the additional leakage 




dproip)Uip)e^^'' 



2\ 



le 



(23) 



which is determined by typical values of U{p). We note that le is of the order of the mean free path, /. The fact that 
le ^ I can be seen from Eqs. ( ^l|)-(p3| ), taking into account that the area A is always bigger than R^. Moreover, for 
short-range disorder, Eqs. (pl|)-(p3D suggest that le ~ I. 
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In order to have t~'^ small, the actual value of U{p) should be suppressed with respect to typical in each box of a 
size k~^. Then the condition > r requires the suppression factor to exceed t/te. The corresponding probabihty 
can be estimated as 

P(r) = (Te/r)^ = {Q/^r^y, (24) 

where N ~ k^A is the number of boxes that "cover" the area of localization of 4'o- In particular, for the ring-shaped 
fluctuations we have N ~ wpok^. We conclude that, due to additional leakage, the prefactor in the density of traps 
is exponentially smafl, i.e., P{t) ~ e:>qi[—k'^A\ii{Q/kl)]. The exponent k^AhiQ should be compared with the main 
exponent, found in Sec. II. Since the above estimate assumed a short-range disorder, kRc <ti 1, the principal exponent 
is 2 {tt^/sY^^ klluQ [Eq. (|l|)]. RecaU that, for short -range disorder we have w ^ k ^{kp^/ InQ)^/^ ^ po ^ m/k, so 
that the product k^A is ~ m^^^ \n~^^'^{Q/kl). The final estimate for the exponent, originating from the prefactor, is 
~ m'*/'^ In^^^ Q. This suggests that traps are the more frequent the smaller is the angular momentum, m. On the other 
hand, in derivation of the main exponent we have assumed that to > InQ. So that the minimal value of (m^ \nQ)^^^ 
is ~ In^ Q. This leads us to the conclusion that the main exponent dominates over the exponent, originating from the 
prefactor. This is because for the main exponent we have klhiQ > In^ Q, since the relation fc^ > InQ was assumed in 
the derivation of the main exponent. 

The above qualitative analysis leading to Eq. ( p^ ) was restricted, for simplicity, to the case of the short-range 
disorder. In the next subsection we present a rigorous calculation of P(r) for the more realistic case of a smooth 
disorder. 

B. Derivation of P(r) 

Denote by p{t) the probability density that the lifetime with respect to additional leakage is equal to r, so that 
p(t) = dP{T)/dT. The rigorous deflnition of this density reads 

p{t)=N j V{U} e-T'^'^UiT-Tu), (25) 

where the normalization constant is defined as 



I nu} e- 



V{U} 



-1 



(26) 



and V{U} is given by 



V{U} = ^JJ dp,dp2U{pi)K{p,,p2)Uip2), (27) 

where K{pi,p2) is related to the correlator K{pi, P2), defined by Eq. (||), as 

J dp'K{p,,p')K{p', P2) =5{pi^ P2) ■ (28) 

Since we are dealing with photons, the value Uq can be expressed through the r.m.s. fluctuation of the dielectric 
constant as Uq = fcpA [see Eq. (||)]. Using the integral representation of the J-function, Eq. ( p5| ) can be rewritten in 
the form 

P{r) = ^fdt e"/- J V{U} e-^'t^>, (29) 
where the auxilary functional, Vt{U}, is defined as 

Vt{U} = V{U} + i^ = JJ dpidp2U{pi)}CtipuP2)U{p2). (30) 
The kernel, ICt{pi, P2), of the functional VtiU} has the form 

ICt{pi,P2) = ^k^^A-^K{pi,p2)+itS{pi,p2) (31) 
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Following the standard procedure of functional integration, we present the fluctuation U{p) as a linear combination 

i7(p) = 5]Ct,^(/.t,^(p) (32) 

where (j)t^fi{p) are the eigenfunctions of the operator fCt with the kernel /Ct 

iCt(j)t.Ap) = j dpiJCtip, pi)0t,^(pi) = At,^(/)t,^(p), (33) 

and At_^ are the corresponding eigenvalues. 

At this point, we note that the operator /Ct is non-hermician. As a consequence, the functions (f>t^f_t form an 
orthogonal basis if the definition of the scalar product is modified to {4>i\4>2) = / dp (j)i{p) 02(p)- To see this, note 
that the kernel ICt{pi, P2) is symmetric with respect to interchange pi ^ p2, ICt{pi, P2) = /Ct(P2, Pi)- Then it follows 
from Eq. ^ 

jj dpidp2 <^t,^i(pi)/Ct(pi,p2) 04,^2 (P2) = 

= J dpi(j)t,p,APl) J rfP2/Ct(pi,p2)(/)t,p2(p2) 

= At,^, J dpi(l)t,t^, (pi) «!>t,M2 (Pi)- (34) 
On the other hand, it also follows from Eq. ( p3[ ) 

jj dpxdp2 (jit,,!^ [pi) /Ct(pi, P2) 0t,^2 (P2) 

= y rfp2 (I)t.,t,2{p2) j rfPl/Ct(P2,Pl)</'t,Mi(Pl) 

= At,tii j dp2(t)t,t,AP2)(t>t,t,2iP2)- (35) 

Comparing the r.h.s. of Eqs. (|3^ ) and (35), we have 

('/'t,^l|0t,M2> = / 0'P'/'t,Ml(P) <^t,M2(P) = '^Ml,A'2- (36) 



Using the expansion Eq. ( |32|), the functional integral Eq. ( p9[ ) reduces to the integration over a/Z complex coefficients 
Ct^fi in the expansion Eq. (|32|). However, due to the restriction that U{p) is real, the pair {ReCt^^, ImCj.^} can be 
"rotated" in such a way that instead of integrals over real and imaginary parts we get a single integral along the real 
axis. 



To proceed further, we note that the real parts of all the eigenvalues At_^ are positive. Indeed, it follows from Eqs. ( |3l| ) 
and (^) that 

„ , _ 1 \ H dpidp2<i)l^[pi)l^{puP2)<i>tAp2) 

- l2fc|A^j //rfp|0..(p)P ■ ^''^ 

Since Av is positively defined, ReAt^^ is positive. The last remark ensures the convergence of all the gaussian integrals 
in Eq. (|^. Thus, we obtain from Eq. (|^ 



1/2 



P{r) = -J rfie-/MM^) . (39) 
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Since the product of the eigenvalues of an operator is equal to its determinant, Eq. (^9|) can be rewritten in the form 

1/2 



det K-o 
det JCt , 



(40) 



It is convenient to present the ratio of determinants in the integrand of Eq. ( pOj ) as a single determinant. This is 
achieved through the following sequence of steps 

det /Co 1 1 

det/Ct det /Co' det /Ct det (^/Cq^'/C, 

1 



det 



1 



det 



1 + 2itkll\'^ KS 



(41) 



where we have used the explicit form ( pl| ) of the operator /Cf. The operators K and S are the integral operators 

with the kernels S{pi, p2) [Eq. ( |22|) ] and K{pi, p2), respectively. We recall that the operator K is the inverse of the 
operator k [Eq. (p8|)]. Upon transformation (^l]), the expression ( ^0| ) for p(t) takes the form 




1 + 2itA:4A2 KS" 



Rc 

° ' 'dot 1 + 2iifc4A2A'S' 



(42) 



To proceed further, we need to analyze the properties of eigenfunctions and eigenvalues of the operators K and S . 



1. Properties of operator K. 

It is easy to see that the eigenfunctions of K are plane waves. Indeed, since the kernel K depends only on the 
difference (pi — P2), we have 

Ke'PP = / dpxKip^ pi)e'PP' = K{p)e'PP, (43) 



so that eigenvalues of K are the Fourier components of the correlator K{p). Thus, these eigenvalues are strongly 
suppressed if p > R^^ ■ For the particular case of gaussian correlator we have 

k{p) = / dpK{p) e'PP = TiRle-P^^'l^. (44) 



2. Properties of operator S. 

Although the eigenfunctions, £,fj,{p), of the operator S are not plane waves, their width in the fc-space is narrow (of 
the order of inverse spatial extent of the function -^o) as it can be seen from Eq. (|2|). To estimate the eigenvalues, 

A^i, of S, defined as 
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J dpiS{p,pMpi) = A^CmIp), (45) 
we first note that these eigenvalues satisfy the following sum rule 

This rule follows from the identity 

^A^= dpdpxS{p,pi)^C^,{p)^ti{Pi) 

= j dpS[p,p) = ^j dp\Mp)?. (47) 
in which the completeness of the set ^^(p) is used, so that 

Y.W^M-S{p-p,). (48) 

The reason why Eq. (^ ) allows to estimate the eigenvalues is their specific distribution. Namely, the first N « k^A 
eigenvalues are almost equal to each other, while the eigenvalues with numbers ^ > N fall off rapidly, faster than 
(iV//i)^. This rapid fall off of A,^ has a simple explanation. The eigenfunctions corresponding to large /x are close to 
plane waves, so that the value of A^^ is determined by the integral of the rapidly oscillating plane wave over the area A 
and is small due to the cancellation. On the contrary, for n < N the eigenfunction changes weakly within the area A. 
Thus, for ^ < N we have Xp, ^ Xq ^ c/{kN). For the particular case of the ring-shaped fluctuations A'' « 2TTpowk'^, 
so that Aq ~ c/{powk^). 



3. Evaluation of the integral (42). 



Short-range disorder. In this case K{p) w K{Q) = 7r_R^, so that the contribution to the integrand of Eq. (|4 



comes from N eigenvalues of the operator S", i.e., det 



when Rr 



0. Then the integration over t in (| 

1 



1 + 2itk'^A'^KS 



1 + 2iTTtk^A'^Xo 



p{t) = 



4t{N/2 
exp 



can be easily performed, yielding 

(27rrAofc(!A2 



N 



, where RcA A 



_ s2~N/2 

2n tXo k^A'^ ) exp 



2 \ kL 



N\ kL 



(49) 



where we have used the large- iV asymptotics of A^! and the fact that Aq ~ c/ {kN). Since the first term in the exponent 
is much larger than the second one, we recover with exponential accuracy the form of P(t), obtained in the qualitative 
consideration. 

Smooth disorder. In the case of a smooth disorder, kRc > 1, it is the fast decay of K{p) [Eq. (|4^)] rather than A^, 
that introduces a "cutoff" of the determinant in Eq. (^). In this case it is convenient to rewrite the determinant in 
Eq. (B3) in the form 



det 



I + 2itk^A'^ k ='[[\^l + 2itXok^A'^Kipf, 



exp 



In ( 1 + 2itXQ k^A^K{pf,] 



(50) 



The sum over ii in the exponent of (|5C|) goes over both projections, Px and Py, of the momentum p. For ring-shaped 
fluctuations it is natural to consider the radial and azimutal components of p. Since w ^ Rc and po ^ Rc the 
contribution to the sum comes only from a single radial component, while the sum over angular component can be 
replaced by an integral. Thus we obtain 
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det 



1 + 2itkl^^KS 



— exp 


Po / 




. Jo 



dp \n{l + 2it\ok'^A^k(j>) 



(51) 



The main contribution to the integral ( pl| ) comes from the domain tAq fcjfA^i^(p) ~ 1, so that p > R^^- It is 

instructive to perform further calculations for more general form of K{p), namely K[p) ^ exp[— (pi?^)"], which 
reduces to Eq. (Q) when n = 2. Substituting this form into Eq. (|l]) and performing two subsequent integration by 
parts, we get 



det 



1 + 2itk',^^KS 



exp 



exp 



Rc \n + l 
ipQ ( n 



Rc \n + 1 



w. 



{n+l)/n 



dw 



(e'"-"'* + i) 



(52) 

^ has a sharp 



where wt = lii{2TrtXoR^kQA ). In the second identity we made use of the fact that the function 

maximum at w ~ Wt- The remaining integral in Eq. ( |5^ ) can be evaluated exactly, J^^dw i^^^_^_^^'2 = — i, so that 
Eq. ( ^ ) takes the form 



det 



1 + 2itk^A^KS 



Substituting this form into Eq. (p^), we obtain 

1 



p(r) 



dt cos{t/T) exp 



exp 



npo 



Po n 



Rc \n + l 



w 



in+l)/n 



2{n + l) 

Evaluating the above integral with an exponential accuracy yields 



^ln("+i)/"(2vrUoi?^fc4A2) 



p{t) exp 



npo J („+!)/„ (2^^XoRX^' 
2(n+l)i?c ^ ' ° 



(53) 



(54) 



(55) 



Since p{t) and P(t) have the same exponential dependence, the final expression for the prefactor, P(t), takes the 
form 



P{t) ~ exp 



npo 



2{n+l)R, 



ln("+i)/" (Q/fc/a) 



(56) 



where r in the argument of the logarithm was replaced by the trapping time u ^Q. For particular case of the gaussian 
correlator m4) the probability that the lifetime with respect to additional leakage is longer than lj~^Q is given by 



P(w"^Q) ~ exp 



(57) 



Note, that for the short-range disorder, [Eq. (^)], the number N was the number of sections with the area ~ /c"^, 
which "cover" the ring-shaped trap. Correspondingly, for the smooth disorder the ratio po/ Rc in the exponent is the 
number of squares with the side ^ Rc that cover the trap. On the other hand, as it is seen from Eq. (p9), the power 
of the logarithm is specific for the gaussian correlator. 



C. Optimal Ring 



Combining the main exponent for the gaussian correlator, Eq. (p!7|), and the corresponding prefactor, Eq. 
obtain 



we 



^UQ)^P{LOcQ)e 



exp 



2,kR, 



In^/^ Q - 34/34-5/3^ 



kl 
JkR^^ 



(58) 
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It is seen that the dependence of the main exponent and the prefactor on m are opposite. The prefactor, reflecting 
the "vulnarabihty" of the ring, faUs off with m, while the main exponent favors large m value for which the Q-factor 
is higher. As a result of the competition between the two tendencies Nm(Q) has a sharp maximum at optimal value 
of m given by 



^opt- 45/4 I p 



In^^/^g (59) 



Substituting m — niopt into Eq. (p8|), we arrive at the final result 

Nopt(Q) = exp 



0-1/2 3/4 jklf''^ J 11/8 Q 



(60) 

If we use a general form of the prefactor Eq. (^6|), then the changes in Eq. (|6^) amount to an additional factor 
2-i3/4(3n)i/4(n+l)-5/4(iin + 2). Also the power of InQ in Eq. ( |60| ) modifies to (5n+ l)/4n. Overall, these changes 
are inessential, so that the result, Eq. (|60|), is rather robust. It shows how the trapping is enhanced due to a smooth 
disorder, when the additional leakage is taken into account. Without the prefactor this enhancement manifested itself 
through the combination (fci?^)^ in the denominator of the main exponent, S',„. With the prefactor, (fci?c)^ is replaced 
by (fci?c)^^^, so that the enhancement is weaker, but insignificantly. 

As a final remark, we note that additional leakage, caused by the scattering out of the plane, can be incorporated 
into the theory in a similar fashion as the in-plane additional leakage. Corresponding changes are outlined in the end 
of Sec. II. Recall also, that for the smooth disorder the suppression of Nm(Q) due to additional leakage is dominated 
by the in-plane scattering processes. 



IV. CONCLUSION 

In the present paper we studied a new type of solutions ]6^j7^ of the wave equation, Eq. (^, in a weakly disordered 
medium. The solutions, dubbed as "almost localized" states, describe a wave which is confined primarily to a small 
ring. In an open sample, of size L much smaller than the two-dimensional localization length 5, the almost localized 
states correspond to sharp resonances, residing in the high-Q ring-shaped cavities, as discussed throughout the paper. 
However, in a closed sample -which would require perfectly reflecting walls- the resonances turn into true eigenstates, 
whose almost entire weight is located at the rings. In this respect the "almost localized" states differ from the 
" prelocalized" states, extensively studied in the context of electronic transport [72-75]. 

We have provided a quantitative theory of the almost localized states and the associated random resonators, and 
pointed out their relevance for the phenomenon of random lasing. We stress, however, that these random resonators 
exist already in the passive medium, and gain is only needed "to make them visible". Moreover, the resonators are 
" self- formed" , in the sense that no sharp features (like Mie scatterers or other "resonant entities") are introduced: 
the model is defined by Eq.(P, which describes a featureless dielectric medium with fluctuating dielectric constant. 



Our study substantiates the intuitive image [^,^ 58 1 of a resonant cavity as a closed- loop trajectory of a light wave 
bouncing between the point-like scatterers. The intuitive picture in [ ^3|j4^ , |5^ assumed that light can propagate along 
a loop of scatterers by simply being scattered from one scatterer to another. Such a picture, however, is unrealistic 
due to the scattering out of the loop. We have demonstrated that the scenario of light traveling along closed loops 
can be remedied. In our picture the "loops", i.e., the random resonators, can be envisaged as rings with dielectric 
constant larger than the average value. The reason why such rings are able to trap the light is that the constituting 
scatterers act as a single entity: only the coherent multiple scattering of light by all the scatterers in the resonator can 
provide trapping. We have also established that correlations in the fluctuating part of the dielectric constant highly 
facilitate trapping. 

We acknowledge the support of the National Science Foundation under Grant No. DMR-0202790 and of the 
Petroleum Research Fund under Grant No. 37890- AC6. 
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